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Existence of Rayleigh Waves Over The 
Heterogeneous Prestressed Medium 

Inder Singh Gupta , Amit Kumar 


Abstract — The frequency equation of Rayleigh waves 
propagating over the free surface of an prestressedorthotropic 
perfectly elastic, heterogeneous semi-infinite medium with 

material properties varying as ^ — ^ , 

P — Pq e a nd principal components of the initial stress 

= S&“.S„ = s 3 V“) ( a > z being the distance 

from free surface) has been derived using matrix method. Then 
for a particular case, for unstressed heterogeneous medium, we 
showed that, there exist single Rayleigh modes which cannot 
propagate below certain cut-off frequency. It is found that in 

case < c < c o w ij ere c denotes phase velocity. ^ is 

Cn 

the constant shear wave velocity of the medium . L is the 
corresponding Rayleigh wave velocity of homogenous medium 

A 

of the same Poisson‘s ratio. On the variation of i.e. 

/l (3 ~-^ - — 

1< < 4 , graph has been drawn between and 1 . 


Index Terms — Rayleigh wave, heterogeneous medium, 
prestressed medium and dispersion equation, phase velocity. 


I. Introduction 

After the poineer work of Lord Rayleigh (1885), the theory of 
surface waves inhomogeneous an isotropic half spaces has 
enjoyed remarkable progress in the 1960’s and 1970’s. Most 
of the works have been summarized in these books of 
Ben-Menahem and Singh (1981), Pilant (1979) and 
Brekhovskikh(1960) for more detailed study. Due to 
complexity of the problem, the study of propagation of 
Rayleigh waves in heterogeneous half-space was not enough. 
It has been studied in heterogeneous medium among others 
Stoneley (1934), Wilson (1942), Newlands (1950), Dutta 
(1963), Hook (1961) and Karlson (1963) . 

Singh (1965) discussed the propagation of Rayleigh type 
waves in an axially symmetric heterogeneous layer lying 
between two half- spaces. The variation of the parameters in 

ft _ & l p 

the layer is assumed to be of the form Ad = 1 ~ az P-d 

l 

_ l where a is constant and " is the distance measured 
from one interface into the layer. He also assumed that the 
vector wave equation for the layer is separable. 

Sidhu(1970) has obtained the frequency equation of Rayleigh 
waves propagating over the free surface of an isotropic, 
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perfectly elastic heterogeneous semi-infinite medium with 
materials properties varying as 

A = A c e az ,ft = pr 0 e az and p = p 0 e az (a > 0), 
He also obtained the solution of the frequency equation in 

i a f ' t /l = ® .• (2) A = ft 

closed form in two cases (1) 

Most of these studies have ignored the initial stresses present 
in the media. Infact, Earth is an initially stressed 
heterogeneous medium .Therefore it is of interest to deal the 
problems on Earth models where initial stresses and 
inhomogeneity factors are taken into account .Biot (1965) has 
shown that an anisotropy is developed due to initial stress 
present in the homogeneous elastic medium. Frequency 
equation of Rayleigh waves in prestressed heterogeneous 
medium is affected not only by initial stresses but by 
inhomogeneity factors also. This study is very useful for 
surface waves analysis. 

Due to complexity of the problem, the study of the 
transmission of Rayleigh waves in heterogeneous half-space 
was not enough. Some researchers have tried to derive the 
dispersion equation for Rayleigh waves in prestressed 
heterogeneous medium e.g. Das et al.(1992), Kakar and 
Kakar(2013) and Abd -Alla et al.(2009).Most of them used 
potential method to solve these problems .Sidhu and Singh 
(1983) and Norries(1983) also have pointed out that potential 
method is not suitable for prestressed media . 

Authors wants to discuss the propagation of Rayleigh waves 
over the free surface of a prestressed heterogeneous 


4= = 


half-space in which 
Po 


and 


g L1 


= € 


az 


and components of normal initial stresses 

-'r r r ? h -^s" 




= e 


and 5 “ 


= e 


) ,( a >0), z 

being the distance from free surface .In the present study we 
investigate the equation of motion with the help of matrix 
method. This method of obtaining the solution of equations of 
motion, can be readily applied to obtain the frequency 
equations of Rayleigh waves in a heterogeneous layer over 
laying in a heterogeneous half- space .These problem are for 
further study. 


II. Basic Equations 

Consider a semi-infinite, perfectly elastic prestressed, 
heterogeneous medium .The materials is either isotropic in 
finite strain or anisotropic with orthotropic symmetry .The 
principal directions of prestress are chosen to coincide with 
the direction elastic symmetry and the co-ordinate axes .The 
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state of prestress is defined by principal components 

b-i-i pSnn Snn r 1 . .... r _ 1 

11 ^ “ oi the initial stress. The general 

equations of motion for prestressed solid in the absence of 

external forces is given by Biot (1965). 

s ij r j ~ — ^ik 0 J jk r j ~ e jk^ ikj —P u i..ttf ^ j 

where ^ is the density, are the displacement 

components and , i indicates differentiation with respect to 

X ■ S - - 

L The incremental stresses are assumed to be linearly 

|0 . 

related to the incremental strains L; through the incremental 

Q-, and Q 2 . 

elastic coefficients J 

5 ll = B ll u x + - w z) > 

S Z2 (^12 P ~) U x + B 22 V y + B 23 W s 

*33 = ( S 12 ~ P ) U x + B 23 V y + B 22 W ^ (2) 

-12 = Q 2 (u y + v x ), 

*13 = Qz ( U z + 

-23 = Ql( W y + V z) ■ 

, S 72 = ^ii ond 

Where, it is assume that 11 

are constants. 

(x,y,z) = {x lr x lr x^, 

{ U,V,W ) = , 

1 


(3) 


du 


etc. 


(4) 


s ij and s ij 


incremental strains. 

Law. 

=^ kk ^ +2 ^iJ ( 5 ) 

Here it is assumed that material properties and initial stress 
components are varying as 

l = JV 1 , ft = P = p"e“ and 5 U = 5 1 " 1 e“,S 22 = S^e cz ,S 33 = S" 3 e “ z , P = 

P ll e° z (a>0) 

( 6 ) 


Then the incremental elastic coefficients becomes 

B u = »ne“ B 22 = Bie az , B u = B\ z e az , B 23 = B? 3 e as 


(7) 


<?i = Q = Qje “ P c = 5 S C S - 5£. 

f B h C h B?. BL , Qi , (??’) S c 5 0 
Where ^ 11 22 12 ’ i3 1 and( llJ oa )are 

elastic coefficients and initial stresses in homogeneous 

U • A A’ J 

orthotropic prestressed medium. L ' are Lame s 
constants and density of material in unstressed state. 

Using (2), (3), (6) and (7) in (1) we get 

d 2 u / „ „ d 2 w ( / n _ P Q \ d 2 v ( n 


e i> = 2 ( u i.j - u ;.d 

_ dlly _ du- 

Ui 'i dxj fUi ^ at dx f 

And P ~ 5ss _ 5l1 . 

Where 

B 1;1 = (2/t -f A)(l + £ 1;1 — 2s z? ') 

B 22 = (2/1 + A) (1- %) 

- . ■. 1 cT -n -n ' j "^11 

^-(1 — ^ll) — ^33 ? 

i 

Ql — U — C P + X)s 22 + - (A — 2 fl) , 

1 1 
Qz =M + -(m + A)(£ii +■ £ 22 )--(A - 2 }i)e 22 . 

J| £■ 

Further ' are Lame’s constants and " L -‘ are the 


are related by Hooke’s 




ftS+ Y)S + a ( <3 ° + ^^ + a l c “ _ “ ,i:i “" a 


2 Jdz 


P Q \ dw (j^ 2 ^ 


2 J dx r dt 2 


| d * w i 




5n 

&y 


Q! ZT - + Cl) S - Bh gf + - P°) 




&U . n nC l _ ^,C e ‘ v '- 

-+ CLB 23 -+a3 22 — -p — . 


( 8 ) 


in. Plane surface waves solution 


Let 


u = u 1 (z)e“ iU,r+: ^- pt) ' 
w = w 1 (z)e _i( * f+3n,_pt) - 


(9) 


where i 

Then displacements in (9) define a plane harmonic wave 
propagating in the direction of the normal to the plane as 
shown in Figure 1 . 



x^ + yti = constant, ( 10 ) 

- ” 2 IT -3 

T T" c = - 

with period ■ , wave length and phase velocity ,CF 
where 

? Jr- 2 n !3 

(11) 

Substituting (9) and (11) in (8) and removing the 
exponential factor , we get 
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dv / 1 

Cuj + A 3 ^v 1 - iaA^Wi -|- i^A 3 - aA 3 — - A t — - = 0, 


A 3 '1% + Dvj -I- iaQ jtiWi 4- ir|A 4 ^ - aQ? ^ = 0, 

du * dv« Sw* 0 z w» 

to {Biz ~ P ) ^i + + Ew t + t£A + irjA^ — td — I = 0, 

( 12 ) 

where 

C = B* ± f 2 + A lV 2 -p a p 2 , 

D = 5 °2 V~ ~ AnO; 2 — p a p 2 , 

E = Q°r +A 2 f 2 -pVA = Qi- F -,A 2 = Q$- F - 


I 3v, 


h,Q . 


T being the column 

Tl T 1 f Tl Tl Tl Tl 

elements lJ 2 ' 3J 4 ' 5 , s and 

6 0 0 0 1 


matrix 


0 0 


0 

0 

c_ 

A 

Aft7 

Qi° 

R° P° 


0 

o 


0 0 
0 0 


A 3 ^rj iaA 2 g -aA 2 


1 

0 

0 


D 

Q\ 


A A 

iarj 0 -a 


0 

1 

AA 

A 

if]A 4 

~qF 


ia(B u n -P")t ia.p 23 vj 


V l 


l 


*A£ iA 4V 

l l 


(15) 


A= 
where 

U 1 = = ^2 ' W 1 = "^7 = ^4 '"^7 = ^5 > = T & 1 


dz 


(16) 

In order to solve (14), we write 

T = p 2Z r 

l, (n) 

=Ae 3Z C 

Which is satisfied for all z if 

[A — 5/]C = 0 

where I denoting six order unit matrix. 


C -r 


(a 2 + 2 (p ± + p 2 ) - 2 V (Pi - P 2 Y ~ IPs ) ' ‘ 


( 20 ) 


i t 

(a 2 +2 (p 1 + p 2 ) + 2 V(P! -P 2 ) Z - 4 Ps) + 1 


A - S 12 + <?; -—,^4 = S? 3 + <?£,£ = B$ 2 . 

(13) 

The set of (12) can be written in the following matrix form 

— =AT. 

dz ( 14 ) 


1. 

s b =~ [V« Z + 4 p 2 - a], 

= ~±[f* + 4 p 2 +a], 

where 

Pi 


2 P 2 9 P 2 2 l{\ _ C s\ 

= ° \1- 2-^ 


with 


Ni 


p D p D 

S = 2T' = ~ and l 7 = 

P Pd Pc- e ii 

( 21 ) 

In case the initial state is unstressed, the medium is 
isotropic and the first order theory of classical elasticity is 
assumed, then it may be shown that (Biot [10]) 

M 2 W 2 

A = A = A + fi Q , C p =a z ,C s = P 2 , p x =tr z -— f p z =tr 2 -—,Pz = a?a z (l-2Y 2 } i 

^11 — ^22 — ^0 ^12 — ^23 — K ?1 “ Ql ~ J^O^l — ^2 “ 

2 ^0- ^ /A]- 2 1*1$ 


Pc 


,P 2 = — and y z = \ - 7-, ) 

Pq ■"4-u 2{I q J 


( 22 ) 
a and 0 


being the constant compressional and shear 
wave velocities in the unstressed homogeneous medium. If 
the medium becomes homogeneous the putting a=0 in (21) 
and (20) becomes. 


The set (18) has a non -vanishing solution vector if and 
only if 

|4 -5/1 - o (19) 

The determinant (19) is a sixth order equation in s and 
gives the following six distinct values, 

t r , if 

[a 2 A 2 (p 1 A p 2 ) - 2-V (Pi - P2) 2 - 4 p 3 ) 

i/ 

(a 2 A 2(p 1 A p z ) A 2y (p 1 -p 2 ) 2 ~ 4p 3 ) ‘ “ - 


s 1 = ^=*J(l-^ 2 =,/p 1 =^l(l-^S S = -aj(l 

s * = 55 = - !) (23) 
The general solution of (14) is, therefore 

f = y e H e s ’A r 

i;=1 1 (24) 

J" 1 (2 

Writing the column vector, ‘ [ in the symbolic 

form as 

^ ^ [ T ^ c i - = (1,2,3, ... . .,6) in (24) and 

equating the corresponding elements ,we get 

Ui = Ti=Z f=i«*‘*Cii 
i\ =T 2 = 2f =1 e SiS C i2 ^w 1 = T 3 = Sf =1 . 

(25) 

We replace s by S * an ^‘ ^ ^ | n Qg) t0 g et 

[>1- Vlfc-l = o 
L 1 JL J} .(26) 

For each i ,(26) gives a set of six homogeneous equations 

. . v C tJ (J = 1,2,. ....,6) 

in six unknowns J . Since the L have 

been obtained from (20), we remark that each of the sets of 


C, , } 

six L -" for any given i must be consistent .We choose a set of 

five from these six E;J fixed, and express them in term of 

sixth in such a way that for none of the values of 1 the 

coefficient determinant vanishes. This happens if we express 

C »0 ~ 2 j , 6 ) . C, = if, (say). 

in terms of ^ E After 

some manipulation, we obtain 
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C„ = 


-k, - , i = 5,6 

77 

t, J , i= % 2,3,4 ' 
0 , i = 5,6 


, (27) 


C,, = 1 im L k 


fl 


, i = 1.2,3, 4 


(28) 


Where 

= ■ _ 123 4 
c (E-Z^-aZsJ 

and 

4 

3 is already defined in (13) (29) 
Substituting for these in (25), we obtain 

u i = 2f=i 

n =fS=, «*<**, 

w i = 


(3“ 3 t4)i^ + Sk = 0. 

Where 

fi = + 5^ 

CO x — — t-rJ-n 3 } fr-hr — — 2 ■' = — ^3 1 .1 

e ~ e zx e yy — e z:z. 

Using (9), (31) in (35), we get 

, (<??+£) *3 + («-£)”**, , K,+ , 

(<?2 + y) 5 + + (<?» - y) + (q° + y) 5 g tf 6 

= 0 , 


- 


, (30) 


[ (<?i -5 3 ° a )tn 3 + (??J a ] 

0 Qi - 0*4 + V - ^=0, (36) 

■(Bi+5^ - (s? 3 5 3 C 3 V > : m 3 5 3l * 3 


[' 


] °+[ 


Appropriate solution to be used from (30) which satisfy the 

i- z -> oo f K t = K, = K- = 0 , 

radiation as are lor 1 *- a , where 

we must ensure that ^ e ( s i^2^ s 3) - ::= - 1 Therefore putting Put ^' B = 


(BSa + ^iW +( fl; C 3 + + ] 

*4 - (4 - Bi° 2 + P a ){ 2 K &=0 


0 in (36) reduce to 


A'l — Kn — Kr — 0 . 

1 b in (30), we get 

= £7 =3 k t e SiZ — K s e** z 

v, =~'1U e v fe* 

1 c-J 1 n 


, (31) 


<?!*! 


[C^ + ^ /2) - 3 + C^.fi C /2) - 3] ^ +[( ?: C + ^ /2) 
o, 

[ C<?1 - ^33 ) ™3 + 

^fii-SS^+Q^ = 0 , (37) 

[ C^i 2 + - (B : c 3 + 5 3 c 3 / _ B? 2 m z s 3 ] K 3 +[ 

( B i2 + - B 2 ° 2 m 4 s 4] tf 4=0 

If we put ^ = <?? 
equation of (37) are identical and hence (37) reduces to 

[W? + f c /2) ^ + W?.« c /2) ™ S] ^ +[( '?? + f c /2) 

*4+0?? 


5 C = 5 I} 

and 11 then first two 


/ 2 ) 4 ] - 0 , 

(Bi + 5^ + (B 2 ° 3+ 5 3 0 3) ^.^ 2 m 353l ^ 


] °+[ 


w, =72?- 

IV. Boundary conditions 
At the free surface boundary conditions are 

Af x - 0 Afy - o, Afz - o at z =0 (32) 

Here we assume that initial stresses are normal and 
S n ^ S X V Z 

-- then ' ^ ' components of the incremental 

boundary forces may be written explicitly (Biot[1961] as 

i/^tju + V-VA + tsii + ^i-^K + fei + ^-VoK r „0 t t V ! D c -CxZ B C W 

kt ( \ , ) ,/ , t c 1 +( B 23 + S 33)V -B 22 m 4 s 4 iC 4 

A /y = (%+V-Vj^ + (%+V,-Vyz>z + (% + Vz-Vd^ ^ J j 

(?2 f ™/2)s 3 f (?2 - i? /2)mg (fl 2 -1 r/Z)s 4 -1 - U u /2)m 4 

( 33 ) (fin + + (B;“ + S? 3 )fJ J - B 2 "™^3 (Su+SuX’ + ^sa +^>j j -BjV 4 s, 

TT 

when are cosines of the angles between the normal and = ^’ (39) 

the jth direction in three -dimensional rectangular Cartesian which is the frequency equation of Rayleigh waves 

propagating over the free surface of heterogeneous 
prestressed medium. 

CASE 1 - when the medium is heterogeneous unstressed 
P c = O') 

(i.e. - then the value of determinate (39) becomes 

Sn + m q s A — m. 


system. Therefore the values of 
if :i = (s 13 + s 33 w y - S n e 3X ) = 0; 

Af y = (s 23 -S 33l u x -S 33 e zy ) = 0; 

Af z = ( S33 +S 33 e-S 33 e aa ) = 0. 


(34) 


The above boundary forces, take the form (using (2)) 

Afy ~ (Cl s 33 ) My+Q 1 ; ? z = 0) at z = o (35) 


s 3 m 3 i — IqCT" 


= 0 


(40) 


A fz = ( B 12 + 

nl, 


CASE 2- when the medium is homogeneous isotropic 

^ Q 

elastic, the value of m from (29) and putting ~ in 

j-'rt 

equation (28) gives the value of 3 and 4 as 
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m ■ — —a 1 1 ~r 
M 


[(<??/% 


(41) 


Using (23) and (41) in (40), we get 


(2- Jr) : -4 1-4 ll-£ = 0 

■ \ or" 


(42) 


giving the frequency equation of Rayleigh waves for the 
homogeneous case as it should. 

Since ^ and 5 we must have The (31) 

reduce to 

»1 _S?=3 K l “* 


[A ,52, 

This equation can have a non zero solution for 
only if 

(QiP* s f (Q*R*„m' 


Kn - k* 


-+‘ IT) + 
and (53) 

.4 — B22 rn s r 

Or 


IT) =0 


=0 

(Q%+P°f 2) 


W 


7 . > L — J E t 


m 5 =0 


(54) 


(43) 


On simplification, the determinately value of (39) becomes equation (50) c >0 and so 


This, however, cannot hold in the present case, since from 

A n „ , > 0 


and 


i i,£) 3 
^4 [ er" 

mE+m^ 3 54 } ]=0, (44) 
where 


N{( V( s a +O + a 2 A c l 


D: „ Isf Is , .. , .. 

L =E- 1 -a \ (1 =3, 4) 

a ^+*V + < b ‘VV 


is real. Hence the only solution of equation (52) is 
+ K ?i + A ' 

=0, giving the tnval solution 
u = v = w = 0 

trom equation (51). 

Here we consider, the values from (22) 

jj 22 _/l c + 2 jU Q _jS.h Z _ A q 

4 = Ci _C? = 


(45) 


We have now to interpret the frequency equation (44) 
subject to condition 


Re( Sl , Sz )>a 
D, D. = 0 


It is clear from equation that the factor ( 
(44) gives only a trivial solution. 
s* 

Put a - (say) 


Sr, S , Sr, 

Substituting the values of ^ and d from(20)in 4 - *=0, 


we get 

[(Pi -VzT - 4 PiF=o 5 

where 


p Zap 2 -:i-2r 2 -' 

c=~= 

■J 

y 

For c to be real, 1-2 >0, i.e. 

^>0 

In this case, from (43), we get 


(48) 


(49) 


(50) 


u i ~ C^3 + 


5 Z 


v i = £«i 

t 

i m 


,(51) 


Using in (38), we get 

Sr, = s, s f , m.n = m* m r 
- and a - 


e , (55) 


a tt' 4 — 
and putting L 

The frequency equation to be studied from (44) becomes 

(r-is* - als -c!sj( : f J s p s P? 1 JCti D +— > 


- + 2 /t c cr-m- 


( 


^0 


3 U 3 * 4 - m{( 

F c 

IhQh + t) s 3 s 4 , (D ls-s,} = 0 

* }{E+ J ‘ rJ 




) 


(46) 

s 3 -0 = 0 of 

(47) 


+2 * }{E+ d 4J (56) 

It is clear from (56), the frequency equation for Rayleigh 
waves depends on initial stress also. 

5 C =5° 

Again we put 11 3 s the frequency (56) becomes 

{ £-k 3 (s a +a)} { E _ is 4 (s 4 + d)}S^ l{ i a + 

,^4^0 “ 5 i C i , . mr cr:5 i C i ( /<o -4 C i) 


2™ ) "" - m[ 


+2 J 


5 C = — o 

Putting 11 ", in above equation we get 

E+V.-^ (57 ) 

In the rest of the analysis we will make use of following 
notations 
a 1 1 ■ 

er= _ L R-V 


'— = ( 

Cs ’ , o'" C li and i/i > 0 f _ 

- sothat =’ v ^ ,">0, (58) 

Cr 

^ being the phase velocity of Rayleigh waves. 
From (20) and (21), using (58), we obtain 

C- - 2 C= p 3 = o 


S n = 


and 


s ■ = 


'(in 1-CT. 


( 59 ) 
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where” 


st = a 2 f 1/1 — 1) 
' 1 + 4 - 


and 


= Q ~ + 4^“ 

Using (59) in (57), we obtained 

= f or (-1 + a) + CL) = 2f 


(60) 


(61) 


From (59), a necessary condition ■ lf >a (here 

s! and si ^ , 

1 - are real) is that □>!. 

From (60) and (61) after some calculation, we get 

a( 4^-z^+i) 


s i=~- 




(62) 


ib ib 

changing (^ -1) to ( - and vise versa in (62) we get 


5n-_ 




(63) 


Since ^ >1, ^ lf in the range 1< v '< 4 

(64) 

Using equation (13), (20), (21), (58) in (45) and again 


putting 


A a — ■ 0 . 


in these equations, on simplification, we get 


” !) 

O, = (65) 

From (60) and (62), we get 


to_ 2 F| 

( 66 ) 


2n 7 ^ ( an J jio/- a 


aB ■ 

tdf — to 00 

it is clear from (66) increase from -- 


as 




3 -f-, 5 


The value of cut off frequency 
aB 


increases from 1 to 4 
below which Rayleigh waves cannot exists is ^ . At v =1 


— I s 




3 -fV 5 


and at 

cj = (3 - yfs)p 2 gives 


Cr ~ Cg. 


D where 


11 e ^ w Rayleigh waves velocity of 

the homogeneous medium of the same Poisson’s ratio. 


V. For numerical calculations 

To derive Rayleigh waves velocity ‘ ’ of the 

homogeneous medium. Put G = ^ in (6), then Lame’s 
constants and density correspond to the homogeneous 

medium. The condition (46) Re ( s iy s z) * a becomes Re 
( ^ 1 ^ 2 ) > ^ , which is possible only for C = 0 iff 

4' 4 ' — 6ip + 1—0 and become indeterminate in 

( 3 +v 5} 


(62), (63)), from here we found the value of 




4 .We 


have drawn the graph between P and on the variation of 
(3-kT; 

from 1 to 4 as shown in figure 2. 
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Figure 2 


VI. Conclusion 

The analysis of the frequency equation carried out for the 

case ^ J shows that, there exist single Rayleigh modes 
which cannot propagate below certain cut-off frequency .It is 
found that speed of Rayleigh waves lies in a limit 
ft < C < c a , w i lere C, fl 


denote the phase and constant 


shear waves velocity in heterogeneous medium and L c is the 
corresponding Rayleigh waves velocity in homogeneous 
medium of the same Poisson’s ratio. 

It is clear from frequency equation of Rayleigh waves in 
heterogeneous prestressed medium depends not only initial 
stress but also inhomogeneity factor (a>0) .Here we discussed 
three cases given below. 

Case 1 . When the medium is prestressed homogeneous the 
frequency equation coincides with the results of Singh and 
Singh (1991). 

Case 2. When the medium is heterogeneous unstressed 

p c = 0 

( ) the frequency equation for Rayleigh waves 

coincide with the results of Sidhu (1970). 

Case 3. When the medium is homogeneous unstressed then 
the frequency equation coincide with the classical results. 
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